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Statement  o-f  Problsm 

This  thesis  derives  and  analyzes  an  approximate  veloci¬ 
ty  solution  that  represents  a  suction  -flow  over  a  flat 
plate.  The  solution  is  derived  from  the  momentum  integral 
equation  in  terms  of  the  suction  and  pressure  gradient  pa¬ 
rameters.  These  parameters  are  adjusted  until  flow  separa¬ 
tion  is  achieved.  Further  analysis  is  done  with  displace¬ 
ment  thickness,  momentum  thickness,  and  friction.  Other 
proven  solutions  are  compared  with  this  analysis. 

Sources  of  Data 

The  information  was  collected  through  a  literature 
review  of  journals,  books,  and  goverment  publications.  This 
included  the  original  studies  dated  from  1945  to  present. 

The  solutions  derived  were  analyzed  through  computer  graph¬ 
ics  and  data  output. 

Conclusions  Reached 

The  approximate  solution  is  in  good  agreement  with  the 
Blasius  solution  for  zero  suction  and  zero  pressure  gra¬ 
dient.  Increasing  suction  velocity  prevents  or  delays  flow 
separation,  and  permits  a  flow  to  withstand  higher  adverse 
pressure  gradient.  The  solution  is  limited  to  a  suction 
parameter  equal  to  -2  and  consequently  an  optimum  suction 
velocity  is  found.  Suction  velocity  decreases  the  boundary 
layer  thickness,  and  adverse  pressure  gradient  increases  the 
boundary  1 ayer  thickness. 
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velocity  pro-file  and  to  observe  the  effects  that  suction  and 
pressure  gradient  have  on  the  approximate  velocity  profile. 

A  dimensionless  velocity  profile  is  solved  in  terms  of 
suction  and  pressure  gradient  parameters.  Pohl hausen 
assumed  a  fourth  order  polynomial  which  is  also  assumed  in 
this  thesis;  however,  different  boundary  conditions  are 
used.  The  boundary  conditions  proposed  by  Torda  are  used  in 
Chapter  3  along  with  the  proposed  fourth  order  polynomial  to 
solve  the  approximate  velocity  profile  used  throughout  this 
thesis.  The  approximate  velocity  profile  is  used  to  analyze 
the  effects  of  suction  and  pressure  gradient  by  observing 
displacement  thickness,  momentum  thickness,  and  shear  stress 


1 


2 


at  the  wal 1 . 

The  -flow  being  analyzed  is  limited  to  a  steady,  two- 
dimensional,  incompressible,  laminar  -flow.  Pressure  gra¬ 
dient  is  limited  to  the  point  o-f  separation.  Another  pur¬ 
pose  o-f  this  thesis  is  to  observe  the  movement  o-f  separation 
pro-file  as  suction  is  applied  to  the  derived  velocity  pro¬ 
file  equation.  The  entire  solution  is  in  dimensionless  form 
where  u  /VJ„  is  restricted  to  remain  between  the  limits  of 
zero  and  one.  If  u  /U  were  to  exceed  the  set  limits,  the 
flow  could  not  be  considered  a  steady  flow. 

In  Chapter  2  the  momentum  equation  is  derived  with  the 
suction  term;  this  equation  then  enables  the  critical  suc¬ 
tion  velocity  to  be  computed.  The  suction  velocity  needed 
to  keep  a  flow  from  separating  is  derived  in  Chaper  4  by  use 
of  the  momentum  equation.  The  derived  equations  from  Chap¬ 
ter  3  (dimensionless  velocity  profile,  displacement  thick¬ 
ness,  momentum  thickness,  and  wall  shear  stress)  are  com¬ 
pared  with  the  solution  of  Blasius  in  Chapter  4  and  analyzed 
for  the  different  suction  and  pressure  gradient  parameters 
which  are  then  tabulated  and  shown  graphically. 


CHAPTER  2 


Derivation  o-f  Boundary  Layer  Equations  with  Suction 

2.  1  Derivation  o-f  the  Boundary  Laver  Equations 
The  Navi er-Stokes  equations  have  been  the  primary  equa¬ 
tions  'for  solving  fluid  flow  problems.  Prandtl  simplified 
the  Navier-StoKes  equations  by  eliminating  some  of  the  terms 
to  develop  the  boundary  layer  equations. 

The  complete  Navi er-Stokes  equations  for  two-dimen¬ 
sional  flow  are: 

du  eu  du  du  ldP  /  a2u  d2u  d2u  \ 

—  +  u—  +  v—  +  w—  * - +  V  I  —z  *  — r  +  — 5-  <  1> 

at  lx  dy  dz  Pdx  dz^  ) 

dv  dv  dv  dv  ldP  /d2v  d2v  d2v\ 

—  +  u —  +v —  ♦  w —  - - ♦  V  I  — =-  +  — =•  +  — —  I  (2) 

dt  lx  dy  dz  Ply  \dx4  ay  dz2/ 

The  boundary  conditions  to  be  satisfied  once  the  boundary 

layer  equation  is  derived  will  be: 

at  y  *  8;  u  *  8 ,  v  »  0 

at  y-  ii  u  =  U* 

Since  the  flow  on  the  flat  plate  is  a  steady  flow  only 
in  the  two  dimensions  of  x  and  y,  the  Navi er-Stokes  equa¬ 
tions  can  be  simplified  because: 


du  av 

>2  2 
a  u  a  v 

du 

av 

1 —  *  w —  ■  0 

;  "  —  m  ■  —  ■  0 

•  — 1  — 1 

f 

=  — 

az  az 

az2  az2 

at 

at 
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The  Navier-Stokes  equations  can  be  reduced  to  boundary 
layer  -form  by  estimating  the  order  o-f  magnitude  o-f  each 
quantity  as  proposed  by  Prandti.  By  letting  the  free  stream 
velocity  UM  be  the  unit  of  magnitude  for  velocity,  the 
terms  of  the  boundary  layer  equations  can  be  put  into  non- 
dimensional  form.*  The  boundary  layer  thickness  &  is 
assumed  to  be  very  small  compared  with  x: 


u  v  p 


x  y  & 

x '  ■  —  j  y'  m  —  ;  i '  *  — 

L  L  L 

where  L  is  the  boundary  length. 

For  the  order  of  magnitude: 

u'  ■  1  since  u  changes  from  6  to  U0 

v'  ■  &  since  y  changes  from  0  to  & 

x'  -  1 

y'  -  * 

Re  -  1  /S2 

Both  u  and  v  are  functions  of  x  and  y. 

As  a  result  of  the  preceding  assumptions,  the  order  of 

magnitude  of  other  terms  can  be  determined  as  follows: 

*u'  a2u' 

—  -  1  j  - 5  -  1 

1  Superscripts  refer  to  endnotes  on  page  78. 
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dU '  1 

d2u'  1 

dy'  i 

dy'2  fc'fc 

dv' 

_  =  c 

d2v'  & 

dx' 

*  2  * 
dx'  1*1 

dv'  & 

d2V' 

dy'  i 

1  5  - 2  " 

dy'z 

i  1 
& 

With  these  newly  derived  terms,  Equations  1  and  2  can  be 
rewritten  into  non-dimensional  -form  as  follows: 

u  ' 


du' 

du' 

x  > •/  -  -  B 

f  d2u' 

d2u' 

dx' 

▼  V  —  m 

dy' 

dx'  Re  ' 

v  dx'2 

dy'2 

<3> 


1  •  1 


For  the  y  direction: 

dv'  Jv '  dP'  1  ( d2v'  d2v'  \ 


<4> 


The  order  of  magnitude  of  the  inertial  terms  equals  the 

order  of  magnitude  of  the  viscous  terms  in  Equations  3  and 

2  2  2 

4.  Since  1  is  smaller  than  1/&  in  Equation  3,  6  u'  /dx' 

2  2 

can  be  neglected.  The  order  of  magnitude  of  1  /Re  is  &  . 

As  a  result  Equation  3  becomes  the  x-direction  of  Prandtl's 
original  boundary-1 ayer  equation: 


du 

du 

ldP  d2u 

c 

1 

♦ 

< 

1 

II 

- ♦  v — ~0 

dX 

dy 

Pdx  by 

Applying  the  same  analysis  of  Equation  3  to  Equation 
4,  shows  that  <1  /P>dP  /dy  —  which  implies  that 

<1  /P>dP  /dy  =  0  since  &  is  very  small.  The  final  forms 
for  the  boundary  layer  equations  are: 


du 

dU 

ldP 

d2u 

u  — 

+  v  — 

s  -  ■ 

— 

+  v - 

(5) 

dx 

dy 

P0x 

dy2 

dp 

— 

=  0 

(6) 

dy 

The 

cont  ii 

nui  ty 

equation  is: 

du 

dv 

—  ♦ 

—  X 

0 

(7) 

dx 

dy 

Condi tion 

s  outside  the  boundary- 1 ayer  are: 

u  * 

Uco  5 

v  * 

0  5 

du  /dy  =  d2 u  /dy2  *  0  5  V  *  0 

Appl 

ying 

these 

conditions  to  the  x-direction  momentum 

equa- 

t  i  on 

.  Equ. 

at  ion 

5, 

leads  to: 

U  <dl)  /dx)  =  -<  1  /p)dP  /dx 
©0  0© 


Now  the  revised  version  of  Equation  5  is: 


du 

du 

dU*  d2u 

u  — 

dx 

*  v  — 
dy 

“ 

u —  +  v  — - 

CD  O 

dx  dy 2 

(8) 

The 

boundary 

1  conditions  are: 

at  y 

«  0 

u  ■  0  ,  v  «  0 

at  y 

■  & 

; 
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The  Prandtl  equations  describe  the  'flow  o-f  fluid  with-* 
out  the  effect  of  pressure  gradient  or  suction,  whereas 
Pohl hausen  solved  the  problem  of  boundary  1 ayer  flow  with 
variable  pressure  gradient,  but  without  suction. 

Pohl hausen's  approximate  method  is  based  on  the  momentum 
equation.  The  momentum  equation  is  derived  from  Prandtl's 
boundary  layer  equation  and  is  often  referred  to  as  the 
von  Karman  momentum- integral  equation. 

2.2  Derivation  of  the  Momentum-Integral  Equation 


Using  Prandtl 's  boundary  layer  equation,  Equation  8, 
and  integrating  it  with  respect  to  y,  from  y  *  0  at  the  wall 
of  the  surface  to  y  ■  •»,  yields: 


n 

du 

duA 

dU 

CD 

♦  - - U_ -  dy  ■ 

V  - 

<9> 

u  \  ** 

ay 

dx/ 

ay 

0 

The  fluid  flow  also  satisfies  the  equation  of  continuity: 

du  dv 
—  ♦  —  -  0 
dx  dy 

By  separating  the  terms  algebraically  and  integrating 
each  term  into  the  continuity  equation,  suction  vq  can  be 
introduced  into  the  equation  as  follows: 
dv  du 

dy  dx 

Integrating  each  term  with  respect  to  y  yields: 


8 


—  dy 
ax 


-  v 


y*0 


Therefore,  the  normal  velocity  component  v<y>  is: 


v(y) 


v< 


u 

0)  -  /  — dy 

Je  >x 


(  10) 


and  v<0)  is  the  velocity  at  the  wait,  or  as  in  this  case, 
v<8>  is  the  suction  term.  Pohlhausen  did  not  apply  suction 
and  consequently  v(0)  =  0. 

In  Pohl hausen's  approximate  solution  the  normal  velocity 
term  used  in  the  momentum-integral  equation  is: 


v<y> 


The  shear  stress  is  known  to  be  T 


<  11) 


-P<au  /dy)  where  Tq  is 


shear  stress  at  the  wall.  To  agree  with  the  momentum  equa¬ 
tion  P  and  P  are  substituted  so  that: 


du 
-pv  — 
ay 


<  12) 


Equation  9  is  rewritten  with  the  substitution  of  the 
new  terms,  Equations  11  and  12.  The  variable,  y,  is  the 
distance  from  the  surface  to  a  point  in  the  free  stream 

3 

i den tified  by  y  ■  h. 


<  13) 
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Integration  by  parts  is  used  to  resolve  the  normal 
velocity  integral.  The  rule  -for  integration  by  parts  is: 


To  de-fine  terms,  let: 


du  du  au 

m  =  /  — dy  ,  dn  ■  — dy  ,  dm  =*  — dy 

] c  ax  ay  ax 

Combining  and  substituting  the  terms  gives: 


Vw  ryw 


au 

—  dy 
IZ  <>x 


ffc 

0  J<d 


au 

u  — dy 
dx 


.  t  —  /  ~dy  ldv 

'0  \ 
at  y  =  0 ,  u  =  0 

at  y  =  h ,  u  =  Ua, 

There-fore  the  right  side  o-f  the  above  equation  becomes: 


*bau 

u®  /  — dy  - 

ax 


au 

u — dy 

ax 


<  14) 


Substituting  Expression  14  into  Equation  13  develops 

4 


the  momentum  equation  into: 


au 


au 


au 


u— dy  -  Uw  /  — dy  ♦  u  — dy  -  U*, - dy  - 


ax 


ax 


'0 


ax 


'0 


dx 


which  can  be  condensed  to: 


au 


au 


dU„ 


2u— dy  -  U„ - U* -  dy  -  -  — 


<  15) 


ax 


dx 


dx 


Since  IL,  »  -f<x)  and  u  “  f<x,y)  ,  and: 


a  dU  du 

-  « —  ♦  Ho¬ 
ax  dx  dx 


le 


and  area  B  are  equal.  Also,  by  definition  Z 2  is  equal  to 
momentum  thickness,  where: 


il 


Figure  2. 1 

Displacement  of  Potential  Flow 


The  momentum  thickness  is  defined  to  be  the  measure  o-f 
the  loss  o-f  momentum  in  the  boundary  layer.  This  is  quite 
similar  to  displacement  thickness  where  potential  -flow  is 
measured. 

The  expression  d/dx,  can  be  taken  outside  the  integral 
in  Equation  14.  The  limit  (0,  «•>  is  independent  of  x  and 
the  general  momentum  equation  becomes: 

d  ■?  dU« 

dx  z  dx  1  P 

The  general  momentum  equation  is  further  expanded  to: 


Pohlhausen  used  this  equation  to  solve  his  flow  problem  with 
variable  pressure  gradient.  In  this  thesis  the  same  problem 


12 


is  solved  but  with  suction.  This  is  done  by  applying  the 
suction  term  to  the  momentum  equation.  As  seen  previously, 
Equation  10  is: 


v<  y)  =  v< 


u 

0)  -  - 
Je  ^ 


dy 


which  is  then  substituted  into  Equation  9  to  solve  for  the 
momentum  equation  with  suction.  It  is  solved  the  same  way 
as  the  momentum  equation  without  suction.  The  final  solu¬ 
tion  is: 

d  ?  dU  * 

dx  *  dx  P 

which  is  expanded  to: 

2  dt_  dll..  T 

U®  “  +  <2i2  ♦  t.  )UW -  ♦  U^vo  -  — 

dx  dx  p 

Therefore  the  momentum  equation  with  suction  is: 

’I 


2  +  —  ♦  — 


<  17) 


The  term  vq  is  the  suction  velocity  <or  blowing) 
through  the  wall.  By  using  the  momentum  equation,  Equation 
11,  the  critical  suction  required  to  Keep  the  flow  from 
separation  is  found.  To  have  separation,  adverse  pressure 
gradient  is  necessary.  Chapter  3  derives  an  approximate 
velocity  profile  that  is  used  in  the  momentum  equation  in 
order  to  find  the  suction  velocity  necessary  to  prevent 
separation.  The  velocity  profile  is  derived  in  terms  of 
suction  and  pressure  gradient. 


CHAPTER  3 


Approximate  Solution  with  Suction 

To  show  the  effects  of  suction  on  a  boundary  layer  a 
method  similar  to  Poh1hausen/s  approximate  solution  is  used. 
However,  according  to  Torda,  Pohlhausen's  boundary  condi¬ 
tions  created  a  problem  when  used  with  his  fourth  order 
polynomial.  A  certain  suction  velocity  value  caused  the 
equation  to  approach  infinity  which  is  not  realistic  in 
actual  flow.  Torda  explained  that  this  problem  would  be 
avoided  by  changing  the  boundary  conditions  and  the  boundary 
layer  equation.^ 

3. 1  Derivation  of  the  Non-dimensional  Velocity  Equation 
Pohlhausen's  non-dimensi onal  velocity  polynomial  is: 

11  2  3  4 

-  ■  a*l  ♦  b^l  +  c^H  +  dV  (  1) 

where  'H  is  a  non-dimensional  term  defined  as  y  /i .  Appen¬ 
dix  A  explains  in  detail  how  Pohl hausen  derived  the  veloc¬ 
ity,  displacement  thickness,  and  momentum  thickness  equa¬ 
tions.  The  boundary  conditions  used  by  Pohlhausen  are  not 
the  same  as  those  proposed  by  Torda.  Pohlhausen's  approxi¬ 
mate  solution  shows  the  effects  of  favorable  and  adverse 
pressure  gradient  on  a  flow.  Adding  suction  to  the  flow 
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changes  the  reaction  of  the  flow.  Suction  from  a  wall 
affects  the  flow  and  separation  for  different  values  of 
adverse  pressure.  The  diverse  changes  of  a  flow  due  to 
pressure  gradient  and  suction  are  shown  in  Chapter  4. 

The  boundary  conditions  used  are: 
at  y  =  8  H  *  0) ,  u  =  0 ,  vq  «  0 

at  y  =  S  (I  *  1)  ,  u  *  U^,  au  /ay  *  a2u  /ay2  =  0 

The  boundary  1 ayer  equation  from  Chapter  2  is: 
du  au  1  Ajp\  /«>2u\ 


u —  ♦  v  — 
ax  ay 


+  v 


<  2) 


P\«Jx/  \  ay‘ 

Outside  the  boundary  layer  u  is  constant  and  equal  to  U0  and 
au  /ax  *  0.  Also,  v  is  zero  causing  Equation  2  to  become: 


10P  dU* 

- -  Uw -  (  3) 

Pax  dx 

Substituting  Equation  3  into  the  boundary  1 ayer  equation, 
Equation  2,  yields: 
au  dU^ 


v0  —  »  IV 
dy 


♦  v 


dx 


<  4) 


The  derivative  of  Equation  4  yields  another  boundary  equa¬ 
tion  at  y  »  0: 


a2u 

“o  2 
ay* 


a3u 


ay* 


The  first  and  second  derivatives  of  Pohlhausen's  polynomial, 
Equation  1,  yield: 


—  *  a  ♦  2bH  +  3cil  ♦  4d*r  (5> 

-J-n 

,2 

U  2 

— r  »  2b  +  6c*\  +  12dil  < 

*rC 

At  y  ■  &  CH  ■  1),  Equations  5  and  6  equal  zero  along  the 
boundary  layer.  Velocity  u  is  a  constant  UM  outside  the 
boundary  layer  'flow  and  is  a  function  of  y  below  the  bound 
ary  layer.  Figure  3.1  shows  the  effects  of  a  velocity  pro 
file  at  the  plate  surface  and  at  the  boundary  layer  thick¬ 
ness.  The  profile  changes  when  suction  is  applied. 


F i gure  3. 1 

Velocity  Flow  Profile 
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Substituting  the  above  boundary  conditions  into  Equations  5 
and  6  respectively  yields: 
iu 

—  =0=a+b+3c+4d  <7> 

ay 

a2u 

— -  -  0  *  2b  ♦  6c  +  12d  ( 8) 


Since  u  *  *t  y  *  l  C*!  *  1),  Equation  1  becomes: 


—  =  l  =  a  +  b  +  c  +  d 


At  the  wall,  where  y  *  0  <*l  ■  0>  ,  the  derivative  of  Equation 
6  becomes: 


-  6c  *  24d-»l 


Applying  the  same  conditions  as  those  used  -for  Equation  10 
to  Equations  5  and  6  yields: 


The  terms  in  Equation  4  and  its  derivative  are  Kept  in  non- 

2  2 

dimensional  form  where  *1  *  y  /i ;  thus  du  /0y,  g  u  /«jy  ,  and 
3  3 

d  u  /}y  can  be  rtwr itten  using  the  chain  rules 
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4U  4u 

4* 

Ju 

*y 

4y 

J1) 

4 

f  0U  1  ^ 

\ . 

Ay2 

4y 

V1* 

/ 

there-foret 

«)2u 

1 

4\ 

ly2- 

‘T2 

1 

& 

1  <>2U  1  <*2U  <)*i 

&  a*Uy  &  <Hl2  «}y 


<  13) 


<  14) 


The  same  procedure  is  used  with  the  derivative  of  Equation  4 
•for  Equations  8,  10,  14,  and  15. 
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! 


v0  r 

-T<2b>  »^r<6c>  (17) 

The  objective  now  is  to  combine  and  rearrange  the 

equations  -formed  so  that  the  coe-f-f icients  of  Equation  1  can 

be  solved.  The  equations  to  be  solved  need  to  be  in  terms 

of  suction  and  pressure  gradient.  According  to  Pohlhausen 

2 

the  dimensionless  parameter,  N  ■  <S.  /v>  dU^  /dx ,  will 

account  for  the  pressure  change  that  is  involved  in  fluid 

flow.  Since  dP  /dx  -  (-PU^JdU^  /dx ,  the  dimensionless 

2  7 

parameter  N  **  <-$  /(MU^JJdP  /dx .  Consequently  U„p  is 
affected  by  the  pressure  change  as  the  flow  proceeds. 

Iglish  theorized  the  development  of  suction  flow  with¬ 
out  pressure  gradient,  and  the  term  he  used  to  account  for 

O 

suction  is  M  =*  vc i  /y.  The  terms  M  and  N  wi  1  1  be  incor¬ 
porated  into  the  newly  derived  equations  of:  u  /U  .  shear, 
displacement  thickness  and  momentum  thickness  To 

solve  for  the  terms  a,  b,  c,  and  d,  multiply  Equation  16  by 
3  /i. 

v0  Uw  dU®  V 

— =  <3a)  -  3 - ♦  -s(6b)  (18) 

i  dx 

Summing  Equations  17  and  18  yields: 
v0  U®  dU®  v 

— 5<3a  +  2b)  -  3 —  -  ♦  — x6<b  ♦  c)  (19) 

l  l  dx 


To  put  Equation  19  in  terms  of  the  coefficients  b  and  c, 
multiply  Equation  9  by  four  and  subtract  Equation  7  from  it. 


This  will  yield 
4  -  c  -  3a  ♦  2b 
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which  is  substituted  into  Equation  19  giving: 

vo  %  *"*>  v 

—z  <  4  -  c)  *  3 - ♦  -r6<  b  +  c) 

i  i  dx  l3 

From  Equation  17,  b  *  3vc  /vQfc .  After  substituting  coeffi¬ 
cient  b  into  the  above  equation,  coefficient  c  is  solved  as: 

3  U^dlt  4ve 

~  *2 
t  dx 

c  *  - - -  <20) 

18k  6v  v0 


Multiplying  the  numerator  and  denominator  of  Equation  20  by 
4  2 

v0£  /v  will  leave  the  coefficient  c  in  terms  of  the  suc¬ 
tion  parameter  M  and  the  pressure  gradient  parameter  N. 


v 
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To  convert  coeff icient  c  into  terms  of  M  and  N,  write: 

-3MM  +  4M2 

c  -  - 5-  <21) 

<  18  ♦  6h  +  hO 

Coe-f  f  i  c  i  en  t  b  is  3vc  /v  ■  3c  /M.  Substituting  Equation 

o 

21  into  the  term  for  b  gives: 

3C-3N  ♦  4M] 

b* - r-  <22> 

<  18  +  <SM  ♦  hT) 

From  the  expression  4  -  c  *  3a  ♦  b,  coefficient  a  is  solved 
as: 


4  -  2b  -  c 


Substituting  Equations  21  and  22  into  the  preceding  equation 
for  coefficient  a  yields: 

24  ♦  <SN  ♦  MN 


<23> 


<  18  +  6M  ♦  M) 

To  obtain  the  coefficient  d,  Equation  7  is  used: 

8  *  a  ♦  2b  +  3c  ♦  4d  <7) 

Again  substituting  Equations  21  and  22  along  with  Equation 
23  into  Equation  7  and  solving  for  the  coefficient  d  yields: 
-6  *  3N  +  2hN  -  6M  -  3M2 

d- - r -  <24> 

<  18  ♦  6M  ♦  n) 

Equation  1  is  now  written  in  its  full  and  complete  form.  To 


help  alleviate  the  congestion  of 


18  ♦  6M  +  M  . 


terms  let  0 
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u  /LT  »  (24^1  ♦  12M*12  ♦  4M2^3  ♦  (6  -  6M  -  3M2)  VS  /  D 
00 

+  (N  ((6  +  M)H  -  9t»2  -  3Mt»3  +  <3  ♦  2M)  V1)  )  /  D 

(25) 

Equation  25  expresses  the  -flow  for  various  suction  and 
pressure  gradient  parameters.  For  various  values  of  M  and 
N,  the  different  separation  locations  are  found.  Whereas, 
Pohl hausen 's  equation,  derived  in  Appendix  A,  shows  flow 
reactions  that  are  dependent  only  on  variable  pressure 
gradients. 


3.2  Derivation  of  the  Shear  Stress  Equation 
The  shear  stress  at  the  wall  where  the  flow  begins  sepa¬ 
ration  is  zero.  Applying  the  same  principle  used  by 
Pohl hausen  and  shown  in  Appendix  A  leads  to: 


To  keep  the  equations  in  non-dimensional  form,  Equation  13 
is  used. 


After  substituting  Equation  23  into  the  above  expressions, 
the  shear  stress  equation  is  derived  in  terms  of  suction  and 
pressure  gradient  as: 

T0*  24  ♦  6N  +  MN 

18  ♦  6M  ♦  M2 


(26) 
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3.3  Derivation  of  the  Pi sdI ac amen t  Thickness  Equation 
Solving  for  displacement  thickness  & ,  is  a  tedious 
process  because  o-f  the  inherent  multiple  computations  of 
approximate  methods.  By  definition,  we  know  that: 


Keeping  non-dimensional  and  using  t|  ■  y  /i ,  gives: 

^  d^l  where  d1'!  has  the  limits  0  to  1. 

Substituting  Equation  1  into  the  above  expression  for  /& 
yiel ds: 

<  1  -  <ai»  ♦  b*l2  ♦  C3  ♦  dVS>d*l 


Integrating  /&  gives: 


h  „ 

(V 

a-n2 

btl3 

cV* 

dH5\ 

1 

(27) 

V 

2 

3 

A 

5  / 

0 

After  expanding  the  above  equation  and  combining  common 
terms.  Equation  27  becomes: 


*1  1 

—  -  — (60  ~  30a  -  20b  ~  15c  -  12d)  <28> 

&  60 

Using  Equation  23  the  term  30a,  yields: 

30a  -  <720  ♦  180N  ♦  30MN)  /  D 


Using  Equation  22  in  the  term  20b,  gives: 

20b  -  <  -  180N  ♦  240M>  /  0 

With  the  aid  of  Equation  21  the  term  15c  yields: 
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15c  =  < -45MN  +  60 M2)  /  D 

and,  the  last  term,  12d,  is  ■found  with  Equation  24. 

12d  =»  <-72  ♦  3 6N  *  24MN  -  72M  -  36M2>  /  D 

Taking  these  four  derived  terms  and  substituting  them  back 
into  Equation  28  yields  the  non-dimensional  displacement 
thickness: 

1  2 

—  =  - <432  -  36N  -  9MN  +  192M  +  36M)  <29> 

i  660 

where  Dis<18+6M+  M2)  . 


<  32) 
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The  second  term  of  Equation  31  is  much  more  involved 
than  the  first  term  since  Equation  25  is  squared.  Equation 
1  could  have  been  assumed  to  be  a  second  or  third  order 
polynomial  to  help  ease  the  numerous  calculations,  however 
accuracy  would  have  been  compromised.  Schlichting  stated 
that  assuming  a  polynomial  with  more  than  four  coefficients 
would  not  increase  the  accuracy  of  the  approximate  solu- 

o 

tion.  It  should  be  noted  that  Schlichting  was  referring  to 
Pohlhausen's  solution  of  flow.  Pohlhausen's  solution  does 
not  involve  suction.  Equation  25  does  include  suction  along 
with  pressure  gradient. 

In  reference  to  Equation  1  and  the  second  term  of 
2 

Equation  31,  <u  /UJ  ,  becomes: 

\2 

—  =»  <a*l>2  +  abH3  +  acV*  *  adH5  ♦  ba3  ♦  b2V*  ♦  bc^l5 


U, 


00 


♦  bdV*  ♦  call4  +  cbH5  ♦  c2il6  +  cdH7  +  daH5 

♦  dbil6  +  dcil7  ♦  d2il8 


Collecting  terms  yields: 

2 

—  I  =  a2il2  ♦  2abil3  +  2acH4  +  2adi»5  +  b2il4  +  2bcil5 


U, 


oo 


+  2bdi>6  +  c2il6  ♦  2cdH7  ♦  d2i)8 


Integrating  the  above  equation  as  described  in  Equation  31 
gives: 
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2bdH7  c2H7  cdH8  d2H9 

+  -  +  -  +  -  ♦  - 

7  7  4  9 

Then  after  applying  the  limits,  the  above  aquation  becomes: 


2bd  c2  cd  d2 

7  7  4  9 


Col  lac  ting  tarms  and  simplifying  them  gives: 

-  <1  / 1240) < 420a2  ♦  430 ab  ♦  504ac 

♦  420ad  +  252b2  ♦  420bc  +  340bd 

♦  180c2  ♦  315cd  ♦  140d2) 

<  33) 

collection  and  simplification  of  terms  Equation 

•  < 1/1240) £  a( 420a  ♦  430b  ♦  504c) 

♦  b( 420c  ♦  252b)  ♦  c(180c) 

♦  d< 148d  ♦  420a  ♦  340b  ♦  315c) 1 


With  further 
33  becomes: 


(34) 
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The  terms  a,  b,  c,  and  d  which  are  Equations  23,  22,  21,  and 
24  respectively  are  then  substituted  into  Equation  34.  This 
task  is  very  involved  and  quite  tedious.  The  work  showing 
the  derivation  of  is  further  shown  in  Appendix  B. 

From  Appendix  B  the  momentum  thickness  equation  is: 

—  =  <58464  -  729N  +  948MN  ♦  40032M  ♦  128 16M2  -  612N2 
& 

-  303MN2  +  6241-1^  -  SSM2!'!2  +  78M3N  +  1872M3 
+  144M4)  /  <  1260(  18  +  6M  +  M2)  2>  <35> 


CHAPTER  4 


Analysis  of  Approximate  Solution  with  Suction 

4. 1  Comparison  of  Approximate  Solution  with  Blasius 

The  previous  two  chapters  went  into  detail  deriving  the 
necessary  equations  to  analyze  suction  flow  over  a  -flat 
plate.  It  was  also  necessary  to  involve  the  effects  of 
pressure  gradient  which  adjusts  the  flow  whether  suction  is 
applied  or  not.  The  equations  in  Chapter  3  had  two  param¬ 
eters  that  were  variable — suction  and  pressure  gradient. 

Some  of  the  equations  are  very  tong  and  involved.  To 
study  the  accuracy  of  this  approximate  solution  for  suction 
flow,  it  is  necessary  to  check  the  approximate  method 
against  proven  data  given  by  Blasius  and  other  authors. 

Blasius"  data  is  based  on  no  pressure  gradient  and  no 
suction.  To  compare  the  thesis  solution  with  the  Blasius 
solution  the  integral  momentum  equation  derived  in  Chapter  2 
is  used! 


Using  the  same  assumptions  as  Blasius,  where  there  is  no 
suction  or  pressure  gradient  to  contend  with,  the  integral 
momentum  equation  is  reduced  to  a  more  simple  form.  For  the 
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problem  considered,  Blasius  assumed  the  free  stream  velocity 
Uw  to  be  constant  outside  the  boundary  layer  and  not  a  func¬ 
tion  of  x;  consequently,  dUM  /dx  equals  zero.  The  revised 
integral  momentum  equation  is: 

di-  /dx  -  T  /pu  2  <  1> 

£  O  ^ 

In  order  to  solve  this  equation  using  the  approximate 
method,  a  velocity  profile  must  be  assumed  which  is  shown  in 
Equation  25  from  Chapter  3.  After  applying  the  conditions 
stated  above,  the  velocity  profile  equation  becomes: 
u  /UM  -  <24*1  ♦  d*l4>  /18  <2> 

The  boundary  conditions  for  this  flow  are: 
y  «  0  (1|  .  0)  ,  u-0  <f  <*l>  «  6> 

y  «  l  (1-  1),  u  •  Uw  <f<n>  «  1> 

where  u  /U„  «  f<*l>  -  <24*1  *  6' I4)  /i8 
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By  definition,  shear  stress  at  the  surface  of  the  plate 


T  is: 
o 


Using  the  chain  rule  and  the  expression,  *1  ■  y  /i ,  Equation 
3  is  revised  as: 


Expanding  XQ  with  the  use  of  the  expression,  u  ■  I4,f  <*l>  , 
yields: 
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Substituting  Equation  2  into  Equation  4  givas  tha  shaar 
strass  at  tha  surfaca  with  no  suction  and  no  prassura  gra- 
diant  appliad.  Equation  5  is  usad  latar  to  dariva  tha 
coeff iciant  of  drag. 


T  «  -HU  <  24  / 18>  <  1  /O 

O  oo 


<5> 


Tha  objactiva  now  is  to  ‘find  tha  valua  of  &  in  Equation  5 
so  i t  can  ba  subs titu tad  into  tha  shaar  strass  aquation. 

Taking  tha  momantum  thicknass  &2  as  spacifiad  in  Equa¬ 
tion  35  front  Chaptar  3  and  applying  to  i  t  tha  conditions  of 
no  suction  and  no  prassura  gradiant,  simp l if ias  tha  momantum 
thicknass  *2  *»* 

*2  /t  -  0.143  <6> 


Substituting  Equations  5  and  6  into  Equation  1  givasi 


— <8.1430 
dx 


1 


PU. 


24  1 


-  mj 


00 


18  & 


Equation  7  is  than  raarrangad  to: 


2  d  24  V\i9 

U*  —  <0. 1430  - - 

dx  18  & 

Simplifying  and  intagrating  Equation  8  yiatds: 

i2  24  Vx 

0.143  - - - 

18  U* 


<7> 


<8> 


(?) 
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Equation  9  finally  yields  the  thesis  boundary  layer  thick¬ 
ness  &  that  can  be  compared  with  the  boundary  layer  thick¬ 
ness  derived  by  Blasius.  The  term  &  is  a  function  of  x, 

< i  ■  Kx>> .  From  Equation  9,  S<x>  is  found  ast 
fc<x>  -  4.318  4  vx  /Um  <19> 

Blasius'  solution  for  boundary  layer  thickness  ist 
tCx)  %  5.6  4  Vx  /Um  (11> 

As  observed  when  comparing  Equations  10  and  11,  the  Blasius 
solution  and  this  thesis  solution  for  boundary  layer  thick¬ 
ness  differ  by  thirteen  percent.  The  non-dimensional  veloc- 
i ty  profile  u  /UM  can  be  solved  by  substituting  Equation  10 
into  Equation  2  where  *1  «  y  /E(x>  .  The  velocity  profile 


then  becomes i 


where  the  limits  are* 

0  <  y  4  /vx  <  5 

Figure  4.1  shows  how  the  approximate  velocity  profile  com¬ 
pares  with  Blasius.  The  comparison  is  close  for  a  no  suc¬ 
tion  and  no  pressure  gradient  flow. 

The  shear  stress  at  the  wall  for  no  suction  and  no 


pressure  gradient  is  found  by  substituting  Equation  10  into 
Equation  5  which  yieldsi 


Figure  4, 1 

Comparison  of  Thesis  and  Blasius  Velocity  Solutions 
No  Suction,  No  Pressure  Gradient 
1 — Thesis 
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318  4  XV  /U* y 

The  above  equation  then  reduces  to: 

XQ  -  -6.389  J~U~7kV  <  12) 

Equation  12  is  the  solution  of  shear  stress  by  the  approxi¬ 
mate  method  with  no  suction  and  no  pressure  gradient.  Tq  has 
a  6.9  percent  error  when  compared  to  Blasius. 

The  displacement  thickness  fcj  is  solved  in  Equation  32 
from  Chapter  3.  Since  is  to  be  compared  with  Blasius, 
again  the  terms  M  and  N  are  assumed  zero.  The  result  is 
Sj  ■  <0 . 4)  S  .  The  term  S  is  equal  to  4.318  4  yx  /U^. 
Substituting  &  into  the  equation  f or  i ^  gives: 

«  1.727  4  VX  /\Jm  (13) 

The  values  o-f  -from  Blasius'  solution  and  Equation  13  are 

almost  exact.  The  same  procedures  used  to  find  X  .  and 

1  o 

u  /Um  are  used  to  find  the  momentum  thickness  Taking 

Equation  39  from  Chapter  3  and  assuming  no  suction  and  no 
pressure  gradient  then  “  <0.143><£>.  Therefore 
becomes: 

i2  -  9 .617  4  V0<  /\im  <  14) 

The  difference  between  Equation  14  and  Blasius'  solution  is 
6.8  percent. 

Other  approximations  have  been  made  and  compared  with 
Blasius,  but  none  of  the  other  solutions  included  suction 
and  pressure  gradient  simultaneously.  Table  4.2  shows  how 
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other  Approximations  made  by  Pohlheusen  end  Prendtl  compered 
with  Blesius  end  the  proposed  epprox imet i on  of  this 
thesis.1®  Figure  4.2  is  e  grephicel  displey  of  Teble  4.1. 
Figure  4.1  comperes  Blesius'  solution  end  the  solution  of 
this  thesis.  The  thesis  solution  is  quite  similer  when  it 
is  reduced  to  e  no  suction  end  no  pressure  gredient 
condi tion . 


B1 esius 

Thesis 

Lineer 

Prendtl 

Pohl heusen 

idx) 

s.e 

4.318 

3.464 

4.64 

5.835 

‘l 

1.72? 

1.727 

1.732 

1.740 

1.7505 

*2 

0.664 

0.617 

0.577 

0.645 

0.687 

*o 

0.332 

0.30? 

0.288 

0.322 

0.343 

f  '<«) 

— 

1.333 

1.0 

1.5 

2.0 

Teble  4. 1 

Comperison  of  the  Integrel  Momentum  Method 
end  the  Blesius  Solution 

Key  for  teble  symbol  si 


&<x) 

"  K1  ‘ 

J  Vx 

‘l 

“  K2  ' 

J  Vx 

l2 

-  k3  . 

J  Vx 

Te  -  K4  (uU.)  J  Ua»  /<xv> 


f  '<«> 


d<u  /U  >  /*i\ 


Figure  4.2 

Comparison  of  Velocity  Approximate  Methods 
with  the  Blasius  Solution 
No  Suction,  No  Pressure  Gradient 

1—  Thesis 

2 —  Linear  Approximation 

3 —  B1 asius 

4 —  Poh 1  hausen 


5 — Prandt 1 
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To  compare  Equation  1?  with  the  Blasius  solution  the 
terms  N  and  M  are  aqua ted  to  zero.  Figures  4.3  and  4.4  show 
the  Blasius  and  thesis  solutions  to  be  equal  at  M  and  N 
equal  to  zero.  Equation  17  will  deviate  ’from  the  Blasius 
solution  as  the  term  his  varied.  This  is  shown  in  Figure 
4.3.  As  suction  is  increased  along  with  Re,  the  coefficient 
of  drag  is  limited  by  the  term  M.  Figure  4.5  shows  that 
when  M  is  approximately  equal  to  -3  the  plot  reaches  its 


maximum  value 


Figure  4.5 

Coefficient  of  Dreg  Versus  Suction 
Parameter  for  N  ■  0 

1 — Thesis 


2 — Edwards 
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As  negative  M  further  incresses  the  calculated  drag 

coefficient  decreases.  For  these  reasons  the  solution 

12 

developed  in  this  thesis  is  not  valid  beyond  M  ■  -2. 

Figure  4.3  shows  how  the  coefficient  of  drag  is  affected  by 

the  suction  term  M  when  pressure  gradient  is  not  present. 

This  thesis  agrees  exactly  with  the  results  presented  by 

13 

J.B.  Edwards  in  an  article  on  laminar  flow.  The  coeffi¬ 
cient  of  drag  between  the  two  studies  begins  diverging  at 
M  «  -2}  therefore  the  value  of  M  in  this  thesis  is  limited 
to  a  value  of  -2  and  greater.  Figures  4.3  and  4.4  also 

confirm  this  fact.  When  v  /U  *  -0.001  at  Re  ■  10*, 

o  ■■ 

which  is  equivalent  to  M  ■  -1,  the  two  plots  agree.  When 

Re  *  10^  and  v  /U  *  -0.001,  M  becomes  -0.316.  This 
o  ** 

comparison  is  quite  reasonable  between  the  graphs  only  up  to 
Reynolds  number  equal  to  5  X  10*  and  for  M  equivalent  to  or 
greater  than  -2. 

Suction  varies  with  Reynolds  number  in  order  to  main¬ 
tain  a  constant  suction  parameter  M.  It  is  observed  that 
too  much  suction  can  be  disruptive  to  a  laminar  flow.  An 

O 

optimum  suction  to  K sap  4  flow  laminar  for  R*  *  10  is  found 

.4 

to  b#  v  /IT  *  1.3  X  10  •  In  summary  only  a  small  amount 

o  a 

of  suction  is  necessary  to  Keep  a  flow  laminar  that  is  nor¬ 
mally  considered  in  the  non-laminar  range  <Re  >  2.5  X  10*) 
The  necessary  vq  /U^  that  will  Keep  a  flow  laminar  is 
effective  as  long  as  i t  does  not  exceed  M  ■  -2. 
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The  solution  deviates  at  values  of  M  <  -2  and  does  not 
reflect  a  true  reaction  beyond  the  limit.  Edwards"  equation 
was  also  limited  at  M  *  -2,  but  he  was  able  to  overcame  that 
limitation  with  another  function. 

4.2  Suction  and  Pressure  Gradient  at  Separation 
From  Chapter  3  the  velocity  profile  was  in  the  form  oft 
u(H>  -  UXaH  ♦  b*l2  ♦  c*t3  ♦  dH4>  ( 18) 

VO 

On  the  right  hand  side  of  Equation  18  each  coefficient  has 

2 

the  term,  D-18*4M  +  M,in  the  denominator  which  is 
shown  in  Equation  25  from  Chapter  3.  The  denominator  has  no 
real  roots)  the  imaginary  roots  are  3  ±  3J  -1.  Conse¬ 
quently  Equation  25  will  not  be  infinite  for  any  value  of  M 

where  M  is  a  function  of  suction  v  .  Torda  observed  that 

o 

the  velocity  profile  would  approach  infinity  if  certain 
boundary  conditions  were  used.  This  also  caused  the  coeffi¬ 
cient  a  to  be  independent  of  suction  which  meant  that  the 
velocity  and  separation  would  not  be  affected  by  suction. 
Torda  later  suggested  the  new  boundary  conditions  used  in 
Chapter  3.  The  suggested  boundary  conditions  did  correct 

the  problem  of  coefficient  a  approaching  infinity  caused  by 

15 

M  in  the  denominator. 

The  separation  point  is  located  where  the  wall  shear 
stress  is  zero.  Shear  stress  is  0<3u  /d'lXd’l  /* )y>  . 
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Taking  Equation  25  from  Chapter  3  and  substituting  it  into 
the  expression,  du  m  0,  yields: 
i u  24  ♦  6N  ♦  MN 

—  -  0  -  - r -  <  19) 

„  <18  *  4M  *  M*>S 

o 

Expanding  Equation  19  and  solving  for  N  gives: 

N  -  -24  /<6  ♦  M>  <20> 

Equation  20  describes  the  pressure  gradient  as  a  func¬ 
tion  of  suction  where  M  ■  vQ(  /v .  Poht hausen's  solution 
to  a  similar  flow  without  suction  was  N  -  -12.  Torda  also 
initially  derived  the  separation  pressure  gradient  as 
N  *  -12,  but  the  boundary  conditions  caused  the  suction  to 
be  ineffective  on  the  flow  and  separation.  Appendix  C  shows 
the  difference  between  the  solutions  of  Pohl hausen  and  this 
thesis  for  separation  point. 

Table  4.2  lists  the  pressure  gradient  at  separation  for 
each  suction  term  M.  The  term  M  can  also  apply  to  blowing. 
Blowing  is  positive  M  and  suction  is  negative  M. 


Table  4.2 

Separation  Pressure  Gradient  with  Suction 


To  have  separation,  adverse  pressure  gradient  must 
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•xist.  As  suction  incraasas  the  fluid  flow  is  abla  to 
ramain  laminar  with  mora  advarsa  prassura  which  maans  that 
tha  saparation  of  flow  can  ba  avartad  or  dalayad. 

Figura  A. 6  dapicts  Pohlhausan's  valocity  prof i la. 
Figuras  4.7  through  4.9  show  tha  valocity  profila  of  this 
thasis  as  a  function  of  tha  suction  tarm  M  and  various 
valuas  of  prassura  gradiant  until  saparation  occurs. 

4.3  Saparation  Valocity 

As  a  fluid  flows  ovar  a  flat  plata,  drag  will  avantual- 

ly  causa  tha  prassura  gradiant  to  bacoma  advarsa  and  this 

will  causa  saparation  to  occur.  To  halp  pravant  or  dal  ay 

saparation,  suction  is  usad.  It  is  possibla  to  find  tha 

amount  of  suction  valocity  raquirad  to  Kaap  tha  flow  laminar 

and  fraa  from  saparation.  Prandtl  was  abla  to  find  tha 

valocity  of  suction  that  would  pravant  saparation  by  using 

id 

tha  momantum  aquation  solvad  in  Equation  11,  Chaptar  2. 
Prandtl  usad  Pohlhausan's  valocity  profila  which  is  a  func¬ 
tion  of  prassura  gradiant  only.  In  this  thasis  to  find  tha 
af facts  of  suction  valocity  and  prassura  gradiant,  Equation 
25  from  Chaptar  3  is  usad  for  tha  valocity  profila.  Tab) a 
4.2  shows  that  with  tha  suction  paramatar  baing  xaro,  tha 
critical  prassura  gradiant  at  saparation  is  N  ■  -4.  Substi¬ 
tuting  M  ■  0  and  N  ■  -4  into  tha  valocity  profila  aquation 
tha  ravisad  aquation  for  valocity  ist 
u  *  2V2  -  il* 


(21> 
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V 

Figure  4.6 

Effect  of  Pressure  Gradient  on 
Velocity  Distribution 
<Poh1h*usen's  Solution) 


Figure  4 .7 

Effect  of  Pressure  Gradient  on 
Velocity  Distribution 

M  -  e 

Figure  notet  (u  /Klw>  >  1,  which  is  unsteady  flow 
occurs  when  N  >  4.3. 
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V 


Figure  4.8 

Effect  of  Pressure  Gradient  on 
Velocity  Distribution 
M  -  -1 

Figure  notet  <u  >  1, 

occurs  when  N  >  3.2. 


which  is  unsteady  flow, 


C(C 
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V 


Figure  4.9 

Effects  o-f  Pressure  Gradient  on 
Velocity  Distribution 
M  -  -2 

Figure  notei  <u  /U^)  >  1,  which  is  unsteady  flow, 
results  when  N  >  4.9. 
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Equation  21  is  used  in  the  definition  of  displacement 
thickness  givings 


Solving  Equation  22  yields  *j  /t  ■  6.5333.  To  verify  the 
above  calculations,  Equation  29  from  Chapter  3  shows  the 
same  results.  To  find  the  momentum  thickness  Equation 
22  is  used  in  the  definition  of  givings 


Equation  23  yields  *2  “  8.1269.  To  check  the  thesis 
results,  Equation  35  from  Chapter  3  was  compared  and  gave 
the  same  answer. 

As  stated  and  assumed  by  Prandtl  the  momentum  boundary 
layer  thickness  is  assumed  constant  which  means  that 
d<&2>  /dx  ■  8 .  In  a  general  flow  over  an  arbitrarily 
shaped  body,  the  free  stream  velocity  will  vary  <UM<x>>. 

It  is  known  that  XQ  /P  ■  6  at  separation,  therefore  the 
integral  momentum  equation  stated  early  in  this  chapter 
reduces  toi 


Solving  for  the  suction  velocity  yieldsi 
v0  -  -<2J2  ♦  *j>dU„  /dx 


<  24) 
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Pr andtl  using  Poh I  hausen's  separation  profile  showed 


that  suction  can  be  used  to  prevent  separation.  Howarth's 
separation  profile  which  is  more  accurate  showed  at  separa¬ 
tion  that! 

v  -  -1 .697  J  -VdU  /dx  <31> 

o  •• 

Howarth  assumed  &  j,  dU*»  /  dx  to  be  constant 

and  worked  the  momentum  equation  with  zero  suction.*^  The 
thesis  Equation  30  does  agree  with  Howarth's  Equation  31 
quite  closely,  and  it  is  clear  that  this  thesis  solution  for 
a  non-suction  flow  does  agree  with  existing  data. 

This  thesis  shows  that  as  suction  is  increased  separa¬ 
tion  is  delayed.  Therefore  suction  will  delay  the  separa¬ 
tion  and  keep  the  flow  laminar  until  the  critical  stage  of 
separation  is  exceeded.  Table  4.3  summarizes  the  data  com¬ 
puted  for  other  various  suction  values.  As  the  suction 
increases  the  boundary  layer  thickness  increases  up  to  the 
point  of  separation.  The  increase  of  adverse  pressure  grad¬ 
ient  for  increasing  suction  explains  the  increasing  boundary 
layer  thickness  at  separation.  The  suction  velocity  permits 
the  flow  to  remain  laminar  longer  thus  making  the  boundary 
layer  thickness  appear  to  increase. 

4.4  Displacement  Thickness 

In  Chapter  3  the  displacement  thickness  was  derived 
in  terms  of  suction  and  pressure  gradient.  As  described 
before,  (j  is  the  mass  flow  deficit  in  the  boundary  layer. 
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M 

N 

‘l 

<D 

♦  262> 

® 

£Kj  3 

<D 

£K23 

0 

-4.0 

0 . S333 

0.7871 

2.000 

-1.5742 

-1 

-4.8 

0.52 

0.776 

2.  190 

-1.699 

-2 

-6.0 

0.5 

0.756 

2.449 

-1.851 

Tabl a  4.3 

Saparation  of  Flow  Data  f or 
furious  Suction  Vatocitias 

Kay  -for  tabl  a  symbols: 

(D  v0  ■  Hl  ♦  2t2KdU„  /dx 
(2)  S  -  £  K  j  3  4-Vdx  /dUw 
(D  v0  -  £K23  4-VdU„  /dx 

Whsn  suction  is  appliad  ths  di spl acsmant  thicknass  should 
dacraasa  sines  ths  mass  flow  assantially  has  bssn  displacad 
clossr  to  ths  fist  plsts  dus  to  ths  suction  fores.  Ths 
prsvious  ssetion  showsd  how  ths  vslocity  flow  compsrsd  to 
Blssius  whsn  ths  suction  and  prsssurs  gradiant  wsrs  zsrosd 
out.  Equation  29  from  Chaptsr  3  givas  ths  di spl acsmsnt 
thicknass  as: 

Sj  -  <432  -  36N  -  9MM  ♦  192M  ♦  36M2)  /60D 

Whsn  M  and  N  sach  squal  zaro,  i }  ■  (0.4)^.  Substituting 

(  ■  4.318  4  vx  /\Jm  into  ths  ij  squation  yislds: 

■  1.7274  Vx  /U^.  This  solution  is  in  sxcsllsnt 
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agreement  with  61 asius' sol u tion  for  Tha  sama  proca- 

dura  is  usad  for  various  valuas  of  tha  suction  paramatar  M 
and  tha  prassura  gradiant  paramatar  N. 

Tha  limits  for  M  and  N  ara  shown  in  Tabla  4.2.  It 
sarvas  no  purposa  to  go  bayond  tha  saparation  point;  a 
staady  flow  is  assumad  so  that  u  /Um  is  not  graatar  than 
ona  and  not  lass  than  zaro.  Figura  4.10  graphically  shows 
how  tha  displacamant  thicknass  dacraasas  as  tha  suction 
valocity  incraasas.  To  countaract  tha  suction  forca,  pras¬ 
sura  gradiant  is  adjustad  until  tha  point  of  saparation  is 
raachad.  Incraasing  advarsa  prassura  gradiant  will  causa 
tha  displacamant  thicknass  to  incraasa.  Figuras  4.11  and 
4.12  show  how  tha  prassura  gradiant  af facts  tha  mass  flow 
for  varying  suction  paramatars.  Figuras  4.10  through  4.12 
ara  assumad  to  hava  tha  sama  boundary  layar  thicknass  as 
Blasius  so  that  tha  af facts  of  suction  and  prassura  gradiant 
on  a  flow  can  ba  comparad. 

4.5  homantum  Thicknass 

Tha  saparaton  prof i la  of  a  flow  can  also  ba  analyzad 
through  momantum  flux  which  is  axplainad  by  momantum  thick¬ 
nass.  Tha  flow  of  a  fluid  naar  tha  surfaca  whara  tha  suc¬ 
tion  is  appliad  may  hava  insufficiant  momantum  to  continua. 
Tha  rasult  may  ba  aithar  a  stall  ad  flow  or  a  turbulant  flow. 
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0.0  l _ I - 1 - 1 - 1 - 1 

0  2  4  6  8  10 

X 

Figure  4.12 

Displacement  ThicKness  Versus  Length  X  (inches) 
for  Different  Pressure  Gradients 
M  -  -2 


If  the  Reynold*  number  is  high  enough  to  Keep  the  41014 
moving,  which  is  normally  the  case,  the  flow  becomes  turbu¬ 
lent. 

Equation  35  from  Chapter  3  is  the  thesis  solution  for 
momentum  thickness.  It  was  shown  how  closely  this  solution 
compared  with  the  Blasius  solution  when  suction  and  pressure 
gradient  were  zero.  Comparing  the  thesis  solution  with 
Pohl hausen's  solution  showed  the  solutions  to  react  similarly 
when  adverse  pressure  was  applied. 

Figures  4.13  through  4.15  show  the  reactions  as  suction 
and  pressure  gradient  change.  Increasing  pressure  gradient 
increases  the  momentum  thickness  initially,  however  the 
thickness  decreases  at  the  half-way  mark  toward  separation 
pressure  gradient.  The  half-way  mark  is  where  the  flow 
profile  changes  shape  and  direction.  The  momentum  thick¬ 
nesses  at  zero  pressure  gradient  and  at  separation  pressure 
gradient  approach  each  other  as  suction  velocity  increases. 
This  indicates  that  extreme  suction  (M  <  -2)  can  cause  a 
flow  to  transition  into  a  turbulent  flow.  How  suction  and 
Reynolds  number  affect  each  other  is  explained  in  Section 
4.1. 

Figures  4.3  and  4.5  show  that  increasing  suction 


increases  the  skin  friction  which  also  relates  to  the  causes 
of  separation. 
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X 


Figure  4.15 

Momentum  Thickness  Versus  Length  X  (inches) 
(Pohlhsusen  Solution) 
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The  values  of  vQ  /Um  and  Reynolds  number  dictate  the  value 
of  M  and  whether  or  not  the  solution  is  in  the  laminar 
region.  Any  value  of  M  less  than  -2  is  considered  in  the 
turbulent  region.  The  thesis  solution  pertains  only  to 
laminar  flow.  Figures  4.13  through  4.15  show  that  momentum 
thickness  is  affected  by  suction  flow,  and  that  at  the 
separation  point  the  momentum  thickness  deviates  very 
slightly. 


CHAPTER  S 


Conclusions 

The  purpose  of  this  thesis  w»s  to  study  suction  flow 
over  *  fist  plat*  snd  to  observe  the  flow  reaction  with 
pressure  gradient.  Ths  objective  was  to  derive  an  approxi¬ 
mate  velocity  profit*  in  t*rms  of  pr*ssur*  gradient  and 
suction  param*t*rs.  This  was  accomplished  by  using 
Poh I  hausen's  assumed  fourth  order  polynomial  with  different 
boundary  layer  conditions  later  proposed  by  Torda.  The 
approximate  velocity  profile  derived  in  Chapter  3  was  used 
to  further  analyze  the  derived  solution. 

The  approximate  velocity  profile  clearly  showed  that 
suction  velocity  is  benefical  in  preventing  or  delaying  flow 
separation.  This  thesis  solution  only  applied  to  small 

0 

values  of  suction  <v_  /U„  *  0.002  or  less)  and  to  Re  ■  10 
or  less.  Figures  4.7  through  4.9  amplified  the  point  of  how 
suction  and  pressure  gradient  can  manipulate  a  flow  reaction 
up  to  the  separation  point.  Xn  observing  the  graphs  for 
drag  caused  by  friction  the  thesis  was  less  accurate  for  the 
value  of  M  less  than  -2  which  is  in  the  realm  of  turbulent 
flow.  The  thesis  did  not  address  turbulence  nor  was  it 
applicable  for  turbulent  flow. 


A*  the  suction  parameter  decreased  (mort  suction  added) 
the  coeff  icient  of  drag  increased.  The  combi  nod  rose  t  ions 
of  suction  and  friction  helped  to  delay  flow  separation. 
Suction  can  fore*  the  flow  to  separate  simply  because  the 
momentum  of  the  flow  exceeds  the  viscous  fore*.  It  was 
observed  that  the  laminar  range  for  Reynolds  number  extended 
as  suction  velocity  was  added.  Normally  the  limit  for  zero 
suction,  laminar  flow  is  Re  *  2.5  X  10*.  Adding  suction 
where  vQ  /UM  »  0.001  would  keep  the  flow  laminar  for  Reynolds 
number  up  to  5  X  10*.  As  the  Reynolds  number  increased  to 
10®,  less  suction,  vQ  /Um  m  1.3  X  10~4  ,  was  needed  to 
keep  the  flow  laminar.  Adding  more  suction  at  this  Reynolds 
number  would  cause  a  non-laminar  situation.  This  shows  that 
suction  does  help  keep  a  flow  laminar  for  increasing 
Reynolds  number  but  with  limitations.  This  set  the  limit  of 
the  approximate  solution  at  M  *  -2.  The  figures  showing 
velocity  profile  clearly  showed  that  at  M  ■  -2  the  suction 
flow  patterns  changed  direction  which  also  indicates  the 
same  limit  to  this  approximate  solution. 

The  solutions  behaved  with  the  predicted  patterns  for 
suction  flow  and  changing  pressure  gradient  within  the  set 
limits.  For  comparing  the  thesis  solution  with  other  proven 
solutions  a  common  denominator  of  zero  suction  was  studied. 
Pohlhausen  showed  separation  to  occur  at  N  «  -12,  the  thesis 
at  N  ■  -4.  The  difference  is  caused  by  the  use  of  different 
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boundary  conditions.  When  Poh I  hausen's  conditions  were  used 
the  solutions  agreed  at  N  ■  **12,  however  the  coefficient 
used  to  find  separation  values  eventually  equaled  infinity. 
Consequently  different  boundary  conditions  were  used  to 
overcome  this  as  proposed  by  Tor da.  The  thesis  agreed  with 
Howarth/s  prediction  for  the  necessary  suction  velocity 
needed  to  maintain  laminar  flow  at  separation  (N  *  -4).  As 
suction  was  added  the  separation  was  moved  further  down¬ 
stream  by  simply  permitting  more  adverse  pressure  gradient 
to  react  on  the  flow. 

The  reason  for  using  a  fourth  order  polynomial  to 
approximate  the  velocity  profile  was  that  it  would  yield  a 
zero  second-derivative  at  the  plate  surface  which  was  a 
boundary  condition.  The  approximate  solution  implies  that  a 
definite  boundary  layer  exists)  whereas  an  exact  solution 
assumes  that  the  boundary  layer  extends  indefinitey. 

The  thesis  solution  has  an  advantage  in  that  a  boundary 
layer  can  be  assumed  and  the  necessary  suction  and  pressure 
gradient  parameters  can  be  substituted  into  the  equation  to 
analyze  the  flow  and  thickness.  Another  acfc'antage  of  the 
thesis  is  that  the  assumed  boundary  conditions  prevented  the 
velocity  profile  from  going  infinite  for  any  suction  parame¬ 
ter.  The  boundary  conditions  proposed  by  Pohlhausen  caused 
the  velocity  flow  to  become  infinite  for  M  ■  -6  which  is  not 
a  realistic  solution,  thereby  Justifying  a  new  set  of 
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boundary  conditions. 

A  problem  encountered  was  the  definitions  of  terms  from 
other  studies.  A  noted  difference  existed  in  the  use  of  the 
terms  *1  and  M.  In  this  thesis  H  ■  y  /i  j  in  other  studies 
H  -  y  /<2S>  ,  ij  -  y  /<  4  2  &)  ,  or  H  -  2y  /t .  Another 
difference  was  with  the  value  for  M{  the  thesis  value  was 
M  •  (v0  /U  >  4  Re,  whereas  in  other  studies 
li  -  (v0  /U  )  4  Re  /2.  Another  noted  problem  was  that  some 
exact  studies  had  their  graphs  and  charts  based  on  parame¬ 
ters  that  were  most  convenient  to  their  studies,  and  it  was 
not  possible  to  duplicate  the  same  parameters  due  to  the 
complexity  of  terms  in  the  approximate  solutions.  It  was 
more  convenient  to  compare  the  thesis  approximate  solution 
with  other  approximate  solutions  due  to  the  similarity  of 
their  methods. 

The  following  conclusions  are  made  about  the  thesis 
approximate  solution  for  suction  flowi 

1.  The  approximate  solution  is  in  good  agreement  with 
Blasius  solution  when  suction  and  pressure  gradient  are 
zero. 

2.  At  separation  the  approximate  solution  agreed  with 
Howarth's  prediction  for  a  no  suction  flow. 

3.  Increasing  suction  velocity  does  help  prevent  or 
delay  flow  separation.  Suction  enables  the  flow  to  with¬ 
stand  higher  adverse  pressure  further  downstream  where  the 


64 

separation  'flow  pattern  is  formed.  More  suction  permits 
more  adverse  pressure  gradient  to  react  on  the  flow. 

4.  The  suction  velocity  required  to  keep  a  flow  lami¬ 
nar  at  separation  point  increased  with  increasing  suction 
velocity.  The  cause  for  this  is  that  as  suction  increases 
the  separation  adverse  pressure  gradient  increases. 

5.  Increasing  suction  parameter  made  the  velocity 
profile  fuller  but  not  as  full  as  Pohlhausen/s  and  other 
solutions.  The  velocity  profile  reached  its  maximum  full¬ 
ness  at  M  ■  -2  before  reversing  the  reaction  of  suction, 
thereby  making  the  approximation  solution  limited. 

6.  Momentum  flux  shown  by  the  momentum  thickness 
varied  with  adverse  pressure  by  initially  increasing  and 
then  decreasing  as  separation  flow  was  reached.  Thesis 
approximate  solution  agreed  with  Pohl hausen's  solution.  In¬ 
creasing  suction  parameter  decreases  momentum  thickness. 

7.  A  possible  solution  to  overcame  the  cancellation  of 
a  term  in  the  assumed  velocity  profile  equation  at  M  »  9  may 
be  obtained  by  assuming  a  fifth  order  polynomial  instead  of 
a  fourth  order  polynomial.  This  may  also  permi t  the  assumed 
profile  to  be  more  full. 


APPENDIX  A 


Flow  Distribution  with  Pressure  Gradient 


Pohl hausen 7 s  approximate  method  was  used  to  solve  the 

velocity  distribution  and  was  based  on  von  Karman's  momentum 

equation.  A  'fourth  order  polynomial  was  assumed  by 

Pohl  hausen  to  represent  the  ‘flows  he  worked  the  equation 

in  non-dimensional  form.  Pohl hausen's  velocity  profile  isi 
2  3  4 

u  /U*  ■  a^l  ♦  b^  ♦  c1!  ♦  d*l  where  *1  ■  y  /t 

The  momentum  equation  from  Prandtl  for  the  x-di recti on  is: 


au 

au 

duw  a2u 

u—  ♦ 

v —  ■ 

—  +  y— i 

( i> 

ax 

ay 

dx  dy 

Boundary  conditions  for  a  steady,  two-dimensional  flow  over 
a  wal 1  are: 

at  y  ■  8  (1  ■  0)  u  ■  0,  v  ■  8, 

at  y  ■  t  <t  ■  1),  u  *  U„,  au  /ay  ■  a2  u  / ay 2  ■  8 
theref ore: 


2  3 

—  -  a  ♦  2bH  ♦  3c  ♦  4d*l 

ay 


<2> 


«)2u 

— 5  -  2b  ♦  6cH  *  12dV 

ay 

At  y 


<3> 


&  which  implies  ’l  ■  1,  Equations  2  and  3  became: 
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au 

—  ■  0  ■  a  ♦  2b  ♦  3c  ♦  4d 

ay 

a2u 

— -  -  e  -  2b  ♦  6c  *  12d 

ay 

2  2 

At  y  ■  0  (I  ■  8) ,  Equation  3  reduces  to  a  u/dy  m  2b. 

Prom  Equation  1,  Uw  < dU^  /dx)  ■  i  /P<dP  /dx>  and 

dP  /dx  ■  0.  The  term,  UM(dUM  /dx)  ,  equals  zero.  Conse- 
2  2 

quently  a  u  /dx  ■  0.  Substituting  the  boundary  condi- 

2  2 

tions  into  Equations  1  and  3,  Knowing  that  a  u  /dy  «  2b, 
coefficient  b  becomes  zero.  Solving  Equations  2  and  3 
simultaneously  yields  that  a  ■  2d.  At  y  «  i  H  ■  1) ,  the 
equation  for  u  /UM  becomes! 
u 

—  -l«a*b  +  c*d 

U«, 

Since  b  ■  0,  the  equation  is  revised  toi 

1  -  a  ♦  c  ♦  d  <4> 

The  equation,  a  *  2d,  now  is  substituted  into  Equation  4 
giving:  1  ■  3d  ♦  c  <5> 

Solving  Equations  4  and  5  simultaneously  yields  d  *  i; 
substituting  coefficient  d  into  the  equation  a  ■  2d,  deter¬ 
mines  that  a  ■  2.  Substituting  the  coefficients  a  and  d 
into  Equation  4,  solves  c  ■  -2.  Therefore  the  coefficients 
are  defined  ast 
a  ■  2 ,  b  ■  0 ,  c  ■  -2, 


d  -  ! 


6? 


It  is  known  that  separation  will  occur  with  adverse 
pressure.  The  condition  for  separation  is: 

4u 

—  -  e 

<)y  y-« 

To  analyze  the  velocity  profile  in  terms  of  pressure  gra¬ 
dient  each  of  the  coefficients  <a,  b,  c,  d>  can  be  written 
in  terms  of  pressure  gradient  p.  Redefining  the  coeffi¬ 
cients  to  a',  b/ ,  c' ,  and  d'  where  a'  ■  a  ♦  p,  b'  ■  b  ♦  pt 
c'  *  c  +  p,  and  d'  *  d  ♦  p;  will  enable  the  velocity  profile 
equation  to  be  solved  in  terms  of  the  dimensionless  coeffi¬ 
cient  p.  The  coefficient  p  is  later  solved  in  terms  of  the 

2 

pressure  gradient  N  where  N  «  /v)dU^  /dx.  N  multi¬ 
plied  by  a  constant  is  equal  to  p.  Substituting  a'  ■  a  ♦  p 
into  Equation  2  *t  y  ■  0  <1  ■  0)  will  initiate  the  steps 
to  solving  the  coeficients  (a,  b,  c ,  d)  in  terms  of  pressure 


gradient.  Due  to  separation  at  the  surface,  a'  «  0.  Since 
a  «  2,  then  p  ■  -2.  At  y  ■  S,  Equation  2  becomes: 

•  ■  2b'  ♦  3c'  ♦  4d'  <6> 

and  Equation  3  becomes: 

8  -  2b'  ♦  6c'  *  12d'  ( 7> 

Solving  Equations  6  and  7  simultaneously,  gives: 

8d'  ♦  3c'  -  8  <8> 

At  y  ■  t  (t  ■  I)  i 
u 

—  ■  I  -  b'  ♦  c'  ♦  d'  <?> 
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Solving  Equations  7  and  9  together  leaves: 

-2  -  Ac'  *  10d'  < 10) 

When  Equations  8  and  10  are  solved  simultaneously,  d'  “  3. 
Since  d'  «  d  ♦  p  and  d  ■  1,  then  p  «  2. 

Using  Equation  10  and  the  term  d'  *  3,  then: 
c '  -  -8 

Since  c/  ■  c  ♦  p,  which  is  rewritten  as  -8  ■  -2  ♦  p,  then: 

p  *  -6 

From  Equation  9,  c'  ■  -8,  and  d'  «  3j  therefore: 
b  '•  6 

Since  b'  ■  b  ♦  p  and  b  ■  0,  then: 

p  *  6 

At  a  certain  adverse  pressure  gradient,  flow  will 
separate.  To  find  the  pressure  gradient,  Equation  2  is  set 
to  equal  zero.  This  occurs  at  the  point  when  shear  stress, 
which  is  proportional  to  ^u  /dy,  equals  zero. 

Equation  2  at  y  ■  0  <*l  ■  0)  yields: 

*u 

—  »  0  m  a" 

Since  p  ■  -2  and  a'  ■  0  •  a  ♦  p,  then: 

0  ■  2  ♦  <-2> 

By  letting  N  equal  the  adverse  pressure  gradient,  separation 
can  be  found. 

N 

a  *  2  ♦  — 

6 
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Therefore,  N  ■  *12  is  the  amount  of  adverse  pressure  that 
will  start  separation  in  a  flow.  The  other  N-dependent 
terms  that  are  associated  with  coefficients  b,  c,  and  d  are: 


N  N  N 

b  - - }  c  -  -2  ♦  —  j  d  -  1 - 

2  2  6 

Substituting  the  above  terms  for  coefficents  a,  b,  c,  and  d 

into  the  equation  u  /UM  derives  Pohl hausen's  approximate 

solution  as : 


'a  4  ”  o  a  a 

—  *  2*1  -  2*1  ♦  *1  ♦  —  <*l  -  3*1  ♦  3*1  *  V)  <11> 

Uoo  * 

r® 

By  definition  the  equation,  <1  -  u  /Um>dy,  is 

the  displacement  thickness.  Since  *1  ■  y  /i ,  then  dy  *  Sd*l. 
Also,  u  ■  when  y  *  i  < “*1  »  l>  .  The  displacement 


thickness  is  rewritten  as: 


(  12) 


Substituting  and  integrating  Equation  11  into  Equation  12 


gives: 


&1  3  N 

&  10  126 

r*> 

By  definition,  S2  uJe  u  1  “  u  /U9)dy.  Using  a  similar 
substitution  for  the  term  dy  as  was  used  for  yields: 


APPENDIX  B 


Solution  of  the  Momentum  Thickness  Equation 


In  finding  the  value  of  the  term  u  /U^  is  squared 

and  integrated.  From  Chapter  3,  Equation  33,  the  relation 


isi 


♦  b<420c 

♦  d(146d 


430b  +  504c) 

252b)  ♦  c  C 180c) 

420a  *340b  ♦  315c)]  Cl) 


The  terms  a,  b,  c,  and  d  have  been  derived  in  Chapter  3 
and  are  defined  by  Equations  23,  22,  21,  and  24  respective 
ly.  To  help  organize  the  mathematical  calculations,  Equa¬ 
tion  1  is  subdivided  into  four  parts. 


Part  One 

Part  One  solves  the  following  expression! 
a<420a  +  430b  ♦  504c)  <2) 

The  terms  of  Part  One  are  expanded  as: 

420a  -  (10080  *  2520N  ♦  420hW>  /D 

438b  -  C-5470N  ♦  7540M)  /D 

504c  -  C-1512MN  ♦  2014M2)  /D 

where  M  -  vei  /v,  N  -  C*1  /y)(dMa»  /dx>  * 
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D  »  18  ♦  ♦  M2 

Substituting  the  expanded  terms  into  Expression  2  gives: 

t  241928  -  7560N  -  2<S208hN  ♦  181440M  ♦  48384M2  *  60 480N 

-  18980N2  -  6552MN2  ♦  45368MN  ♦  12094M2N  ♦  18888MN 

-  3150f*l2  ♦  1892M2N2  ♦  7560M2N  ♦  2816M3N  1  /  D2 

Expression  2  is  -further  reduced  by  combining  similar  terms 
to  give: 

t  241928  -  15120N  ♦  29232MN  +  181448M  ♦  48384M2  -  18988N2 
-  9702W2  ♦  1965^^  -  1092M2N2  +  201«5M3N  1  /  D2  C3> 

Part  Two 

Part  Two  solves  the  -following  expression: 
b<428c  ♦  252b)  <4> 

Expanding  the  terms  of  Expression  4  gives: 

428c  -  C-1260MN  ♦  1488M2)  /0  <5> 

252b  -  (-22A8N  ♦  3824M)  /D  (6) 

Substituting  Equations  5  and  6  into  Expression  4  yields  the 
result  of  Part  Two  which  is: 

t  U340MN2  -  30240M2N  +  204 12N2  -  54432MN  ♦  362886M2 

♦  28160M3  3  /  D2  (7) 


Part  Three 

Part  Three  solves  the  expression: 


c<  1880 


Alter  expending  end  combining  similer  terms,  the  solution 
for  Pert  Three  becomesi 


t  I62wh42  -  4320M3N  ♦  2880M4  ]  /D2  <8> 

Pert  Four 

Pert  Four  is  the  expension  of  the  expression: 
d< 140d  +  420e  ♦  360b  ♦  315c)  (9) 

Teking  eech  term  from  Expression  9  end  expending  it  gives: 
140d  -  (-180  *  420N  ♦  280MN  -  840M  -  420M2)  /D  (10) 

420 e  «  <10080  ♦  2520N  +  420MN)  /D  (11) 

360b  -  < -3240N  ♦  4320M)  /D  (12) 

315c  ■  <— 945MN  ♦  1260M2)  /D  (13) 


Substituting  Equetions  10,  11,  12,  end  13  into  Expression 
yields  : 

I  -55440  ♦  1800N  ♦  1470MN  -  20880M  -  5040H2  ♦  27720N 

-  900N2  -  735hW2  ♦  104401*1  ♦  2520M2N  ♦  18480MN 

-  490M2N2  ♦  6960M2N  ♦  1680M3N  -  55440M  ♦  1800MN 

♦  1470M2N  -  20880M2  -  5040M3  -  27720M2  ♦  900M2N 

-  600MN2  +  735M3N  -  10440M3  -  2520M4  1  /  D2 

After  combining  ell  common  terms,  the  expression  of  Pert 
Four  is: 

t  -55440  ♦  29520N  ♦  32190F-N  -  76320M  -  53640M2  -  900N2 

-  1335hW2  ♦  11850M2N  -  490M2N2  ♦  2415H3N  -  15480M3 

-  2520M4  ]  /  D2  (14) 


l 
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Adding  Expressions  3,  7,  8,  and  14  and  then  multiplying  them 
by  < 1  /1260)  gives  the  expanded  version  of  Equation  It 

d*i  -  [  186480  *  14400N  ♦  6990MN  *  105120M 

♦  31032M2  4  612N2  ♦  303MN2 
4  Sen2!'!2  4  lllM^  4  4680M3 
4  12 64M2H  4  360M4  1  /  12600  <  15) 

By  definition  the  momentum  thickness  equation  ist 


Since  Equation  15  from  above  and  Equation  32  from  Chapter  3 
are  terms  in  the  momentum  thickness  equation,  the  new  equa¬ 
tion  for  /i  with  suction  and  pressure  gradient  ist 

<648  4  36N  4  9MN  4  168M  4  24M2) 

&  60  <18  4  611  4  M  > 

-  C  <186480  4  14400N  4  6990MN  4  105120M  4  31032M2 
4  612N2  4  303MN2  4  1266*1^  4  38mV  4  1 1  lM^ 

4  4688M3  4  360M4)  /  <1260(18  4  6M  4  M2) 2>  ) 

After  combining  similar  terms,  the  non-dimensional  form  for 
the  momentum  thickness  equation  becomes t 
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—  -  (58464  -  729N  ♦  948MN  ♦  48032M  +  128 16M2  -  612N2 
& 

-  303MN2  +  624M^N  -  SSM^2  +  78M3N  +  1872M3 
♦  144M4)  /  <1260(18  +  6M  ♦  M2)  2> 


(  16) 


APPENDIX  C 


Derivation  of  Shear  Stress 


To  find  separation  pressure  gradient,  Poh I  hausen  first 
derived  the  shear  stress  equation  and  set  it  equal  to  zero. 
To  have  separation  the  shear  stress  must  equal  zero.  The 
shear  stress  equation  ist 

(du  \  iu  a*! 

—  rn  H -  (note*  a*»  /dy  «  I  /i> 

«>y  )0  a*i  ay 

Taking  Equation  11  from  Appendix  A  and  substituting  it  into 

the  shear  stress  equation  gives: 

par  34n  234-1 

X  » - u  <2*1  -2*1  ♦*!)♦—  <*1-3*1  ♦  *1  -  *1  ) 

0  i  a*lL  *  6  J  a 


The  above  equation  reduces  to* 


Rearranging  Equation  1  gives* 


V  N 

■  2  ♦  —  <  2) 

PU„  6 

Equation  2  is  Poh 1  hausen's  solution  for  shear  stress  in 
terms  of  pressure  gradient  only.  Equation  26  from  Chapter  3 
shows  the  shear  stress  as* 
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Tl  24  +  4N  +  MN 

- -  - -  <3> 

PU*  18  ♦  4M  ♦  M* 

Equation  3  is  a  function  of  both  pressure  gradient  and 
suction.  Pohl hausen" s  solution  is  a  function  of  only  pres¬ 
sure  gradient.  To  find  the  pressure  gradient  at  separation 
where  the  shear  stress  is  zero,  Equations  2  and  3  are  used. 
These  two  equations  each  equal  zero  at  separation,  therefore 
their  pressure  gradients  at  separation  can  be  found. 
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